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Abstract 

We investigate in detail the associated production process pp A'^Z'^ + X within 

the framework of the minimal supersymmetric standard model (MSSM) at the CERN Large 
Hadron Collider (LHC), considering both contributions from the Drell-Yan and gluon fusion 
subprocesses. We focus on the deviations from the general two-Higgs-doublet model (2HDM) 
arising in the MSSM. We also discuss the contributions of the two J^Z'^ associated production 
subprocesses in the MSSM at the LHC, and analyse the dependences of the total cross section 
on neutral CP-odd Higgs boson mass vtla and tan /3 in the mSUGRA scenario. We find that 
the contribution from loop mediated gluon fusion subprocess can be competitive with that 
from the Drell-Yan subprocess in some parameter space. 
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I Introduction 

The minimal standard model(MSM) ||l| [Q] has been proved by ah precise experimental data that 
the MSM is a very successful model of particle physics. But until now the symmetric breaking 
structure of the electroweak interactions has not yet been directly explored experimentally. So 
the exploration of the SM Higgs boson is a major goal of the present and future colliders. As we 
know, any enlargement of the Higgs sector beyond the single SU{2)l Higgs doublet of the MSM 
necessarily introduces other neutral Higgs bosons and charged Higgs bosons. Like the general 
two-Higgs-doublet model(2HDM), the minimal super symmetric standard model (MSSM) |^ |Q] 
requires the introduction of two Higgs doublets in order to preserve supersymmetry. These 
two Higgs doublets predict some more elementary Higgs bosons: one CP-even neutral Higgs 
boson(//''), one CP-odd neutral Higgs boson(y4*^) and two charged Higgs bosons (H^), which are 
absent in the MSM. Any experimental discovery of these non-SM-like Higgs bosons will be the 
direct verification of these extended versions of the Higgs sector. Therefore, the study of various 
production mechanisms of the non-SM-like Higgs bosons at the present and future colliders is 
well motivated. 

Searching for the non-SM-like Higgs bosons and studying their properties at the future multi- 
TeV hadron colliders, such as the CERN Large Hadron Collider (LHC), are possible as expected 
by supersymmetric (SUSY) theory |^] The gluon fusion mechanism gg — > = , , A^) 
provides the dominant production mechanism of neutral Higgs bosons at the LHC in the entire 
relevant mass range up to about 1 TeV for the small and moderate values of tan (3 in the MSSM 
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1^]. The heavy neutral Higgs boson can be also produced in pair {A^A^,A^h^,A^H^) at the 
LHC, if it is kinematically allowed ||^]. Studying the process of a heavy charged Higgs boson 
associated with W boson is another attractive way in searching for the bosons, because the 
Vl^^-boson's leptonic decay may be used as a spectacular trigger. The calculations of the heavy 
production associated with boson at a future electron-positron collider can be found 
in Refs. § [^] |10| |11] The complete calculations of the H^W^ associated production 

at hadron colhders both in the 2HDM and the MSSM are given in Refs. IH |l5| [jl6| 0. 
Analogously, A^Z^ associated production would also be an efficient way in searching for the 
heavy neutral CP-odd Higgs boson A^. Although the A^ boson can be produced in pair at future 
colliders |0] fl^ , the A^Z^ associated production will be the kinematically favored mechanism 
to produce A^ Higgs boson for the heavy A^ Higgs boson. And again the leptonic decay of Z^ 
maybe benefit for triggering the Z^A^ associated production events. The calculations of the 
Z^A^ associated production at a electron-positron collider were presented in Refs. pO], at 



a muon collider in Ref. [21| and at a photon collider in Ref. [22|, respectively. And Chung 
Kao gave the calculation of A^Z^ associated production via gg fusion including only quark loop 
diagrams at the SSC p^ ]. 

In this paper we concentrate on studying the A^Z^ associated production at the LHC in 
the MSSM, considering both subprocesses qq A^Z^ and gg A^Z^. In the calculation of 
the loop mediated process pp ^ gg ^ we compare and discuss the cross sections in 

the general two-Higgs-doublet model (2HDM) and the MSSM. In section II, we present the 
calculation of the processes pp ^ qq ^ A^Z^ and pp ^ gg ^ A^Z^. Numerical results and 
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discussion are given in section III. There we use the MSSM parameters constrained within the 
minimal supergravity (mSUGRA) scenario |24]. Finally, a short summary is given. 

II Cross Section Calculation 

In our calculation we use the t'Hooft-Feynman gauge and adopt the dimension regularization 



scheme in the general 2HDM and the dimensional reduction (DR) scheme |25] in the MSSM. 
In the loop diagram calculation we adopted the definitions of one-loop integral functions in 
reference |2^]. The numerical calculation of the vector and tensor loop integral functions can 



be traced back to scalar loop integrals as shown in the reference [ 27 1. The Feynman diagrams 
and the relevant amplitudes are created by FeynArts package automatically [^] . The numerical 
calculation of the loop integrals are implemented by using Mathematica programs. 

II. 1 Calculation of the subprocess qq — > y4°Z° + X 

We denote the A^Z^ associated production via Drell-Yan subprocess as 

qipi)+qiP2)^Z''{h)+A'{k2), (2.1) 

Due to the feature of the Yukawa coupling that the coupling strength between quarks and 
Higgs boson is in proportion to the correspondent quark mass, the cross sections of subprocesses 
qq A^Z^{q = n, d, s, c) should be much smaller than those of the subprocesses tt{bb) A^Z^. 
Considering the fact that the luminosity of top (anti-top) quark is much lower than that of 
bottom (anti-bottom) quark from a proton, we conclude that the cross section of the process 

— > qq ^ A^Z^ + X is approximately equal to the cross section of pp — > bb A^Z^ + X. 
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Therefore, in this paper we consider only the contributions from the pp ^ bb ^ J^Z^ + X 
process. 

The Feynman diagrams of the subprocess 66 A^Z'^ at the lowest order are depicted in 
Fig.l. The differential cross section of the subprocess 66 A^Z^ can be expressed as 

d^bb = dP2f ^ Yl ^) + (2.2) 

spin 

where the summation is taken over the spins of the initial and final states, and dP2f denotes the 
two-particle phase space element. The factor 1/12 in above equation comes from the averaging 
over the spins and the colors of the incoming partons. The matrix element ^(<j) represents the 
amplitude of the h^{H'^) exchanging s-channel diagrams(shown in Fig.l(a)), A(i,^ corresponds to 
the amplitude of u- and t-channel diagrams (shown in Fig. 1(b)). The Mandelstam kinematical 
variables are defined as 

s = (Pl+P2)^ i={pi-kif, u={pi-k2f, (2.3) 

By using the relevant Feynman rules, we obtain the explicit expressions of these amplitudes: 

A f -\ (47ra)mfe sina {h + Pi + P2T 

A(a){s,t,u) = - .22 cos /3 - a — — e^(fci -u F~ 

4:sfycfy COS /? s - uif^ + nihT hi 

(47ras)m6 . \'^'^^^ uti. \ \-i \ 1 w (^1+^1+^2)'' ,^ 
J^-^— sin(/3 - -^e^l^i) \vm)u{px)\ — (2.4) 



A^i,)is,t,u) = ^tan/3e'^(/ci)^ 

(47ra)mb ^ „/, n 1 
^ tan ^e^ {h ) ^ 



where and mw represent the masses of bottom quark and W boson, respectively. 



II. 2 Calculation of the subprocess gg — > + X 

We denote the APZ^ associated production process via gluon fusions as 

«) + 5(P2, /?) ^ Z\k^) + A\k2), (2.5) 

where a, /3 are the color indices of initial gluons. As the subprocess gg — Z^A^ is loop-induced, 
the one-loop order calculation can be simply carried out by summing all unrenormalized reducible 
and irreducible one-loop diagrams and the results will be finite and gauge invariant. We denote 
(Tg^^^ and CF^g^^^ as the cross sections in the framework of the general 2HDM and the MSSM, 
respectively. The former is contributed by the Feynman diagrams involving only the quark loop 
diagrams (shown in Fig.2) and the latter involves the contributions of both the quark and squark 
loop diagrams (shown in Fig. 2-3). The possible corresponding Feynman diagrams created by 
exchanging the initial gluons or the two final states, should be also included in Fig.2 and Fig.3 
and involved in our calculation . 

We can see that each Feynman diagram in Fig.2 and Fig.3 contains one interacting vertex 
between (s)quarks and a Higgs boson. Due to the feature of the Yukawa coupling as we mentioned 
above, we can consider only the diagrams which involve the third generation (s)quark in the 
calculation of the subprocess gg J^Z^ . The cross sections of the subprocess gg J^Z^ in 
the general 2HDM and the MSSM can be expressed respectively as 

(I, t, u) + ^jg (s, f, u) • • • + (s, t, u)f (2.6) 
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where the summation is taken over the spins and colors of the initial and final states, and dP2f 

(i) 

denotes the two-particle phase space element. A^-^ represents the amplitude of the diagram 
of Fig.i(j). The factor 1/256 results from the averaging over the spins and the colors of the 
incoming partons. 

II. 3 Cross section of pp — > + X process at the LHC 

With the cross sections of the related subprocesses, the cross section of parent process pp 
A^Z^ + X at the proton-proton collider LHC can be obtained by doing the following integration, 

dCi 



CTij = / 

J (m 



dr^^aijis = Ts) 



(2.7) 



where 



dr 1 + 6, 



^ dx^ 



Xi 



fi/p{i,xi,Q )fj/p{j, —,Q 



+ 



fj/p{j,xi,Q )fi/p{h—,Q ) 



(2.i 



In Eq.( |2.7| ) ^/s and \/I are the colliding proton-proton and parton-parton c.m.s. energies respec- 
tively. The notation aij represents the cross section of the parent process pp ij — > A^Z^ + X. 
dCij/dr is the luminosity of incoming partons where i, j can be b, b and g, t = xi X2- mz and 
iTiA represent the masses of boson and A^ Higgs boson. The definitions of xi and X2 can 



be found in Ref.|29]. In our calculation, we adopt the CTEQ5 parton distribution function 
and take the factorization scale Q to be The Eq.( |2.7| ) can be rewritten as 



I d\/laij{s)Hij{s) 

Jmz+mA 



f2.c 



where 



2dxiv§ 

XlS 



fi/p{i,xi,Q )fj/p{j,—,Q ) 



+ (i ^ j) 



(2.10) 



When ij = gg, crl^^^ represents the cross section of the parent process pp ^ gg ^ J^Z^ + X 
contributed only by quark loop diagrams shown in Fig. 2, while (^gg^^^^ represents the cross 
section contributed by both quark and squark loop diagrams shown in Fig. 2 and Fig. 3. In the 
next section we shall take different input data sets to demonstrate the production rates of the 
parent process pp — > A^Z^ + X. The numerical results of cr^j^^^^ and (^gl^^^^ would show the 
importance of squark loop diagrams. The total cross section of pp A^Z^ + X at proton-proton 
collider should be the summation of fTj,^ and agg. Quantitatively comparing the with agg 
will help us to know in which part of the parameter space the contribution of gluon-gluon fusion 
process is dominant. 

Ill Numerical result and discussion 
III.l Input parameters 

In the numerical calculation, we take the SM parameters as: mt = 174.3 GeV, rrih = 4.2 GeV, 
mz = 91.187 GeV, Tz = 2.49 GeV [^], and take the supersymmetric parameters being con- 



strained within the minimal supergravity (mSUGRA) scenario |24]. In this scenario, only five 
sypersymmetric parameters should be inputed, namely M1/2, Mq, Aq, sign of fj, and tan /3, where 
Mx/2) and Aq are the universal gaugino mass, scalar mass at GUT scale and the trilinear soft 
breaking parameter in the superpotential terms, respectively. In this work, we take Mi/2=120 
GeV, Aq=300 GeV and > 0. Mq is obtained quantitatively from the input niA value. All 
other MSSM parameters are determined in the mSUGRA scenario by using program package 
ISAJET 7.44. In this program, the renormalization group equations (RGE's) l3^ are run from 
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the weak scale mz up to the GUT scale, taking all thresholds into account in order to get the 
low energy scenario from the mSUGRA. It uses two loop RGE's only for the gauge couplings 
and the one-loop RGE's for the other supersymmetric parameters. The GUT scale boundary 
conditions are imposed and the RGE's are run back to mz-, again taking threshold into account. 

Here we give some comments about the choice of the decay width values of CP-even neutral 
Higgs bosons and . We know that some of the Feynman diagrams (shown in Fig. 1-3) have 
s-channel and propagators, which have analytical expressions respectively as 



1 s - im\ - rrihTh 

s - ml + imhTh {s - ml)'^ + mlTl' 



(3.1) 



1 _ s-mjj- imnTH 2) 



s — rn^j + iniH^ H {s — ni^jY + ■m\jT']j^ 
It is clear that the cross sections of the subprocess should related to the decay widths of and 
H^. In this work the input parameter rriA is taken in the range of 200 GeV to 650 GeV. Then 
we have the following constraints in this parameter space, 

rriH ~ w-A, "rrih < 150 GeV, (3.3) 
and by using the package HDECAYpI] in the MSSM, we find 



Th, Th < 10 GeV. (3.4) 

Because \/l > niA + mz, we get {s — mfj)^ ^ m'j^T^. The propagator of boson can be 
expressed approximately as 

s — mfj — imu^H 1 



(3.5) 



It is obvious that the bigger the \fi is, the less sensitive the cross section to the decay widths 
of neutral Higgs bosons and is. Therefore, we choose Th = T/^ = 10 GeV in our 
numerical calculations. Actually, our final numerical result of the cross section of the process 
pp ^ A^Z^ + X at the LHC, shows also that it is not sensitive to the choice of these two decay 
widths. 

III. 2 Discussion and analysis 

The figures in Fig.4, Fig.5 and Fig.6, show the cross sections (or differential cross sections) of the 
process pp ^ gg A'^Z^ + X at the LHC as the functions of the CP-odd Higgs boson A° mass, 
the ratio of the vacuum expectation values tan (3 and the transverse momentum pT, respectively. 
The curves of the cross sections (or differential cross sections) involving the contributions from 
quark loop diagrams (in the general 2HDM) and quark+squark loop diagrams (in the MSSM) 
are depicted separately on these figures for comparison of the cross sections in these two models. 
And in these three figures the full- lines are for cr^^^^ (or dag^^'^ /dpr), the dotted- lines are 
for af/^^ (or daf/^^ /dpr). 

Figure 4 shows the relationship between the cross section of the parent process pp ^ gg ^ 
A^Z'^ + X and niA with the colliding energy ^/s = 14 TeV. The input mSUGRA parameters 
are set to be the typical values mentioned in the last subsection (i.e. Mi/2=120 GeV, Ao=300 
GeV and /x > 0. Mq is obtained quantitatively from the value), and tan/3=2, 7 and 
32, respectively. From this figure, we find that in some parameter space the scalar quark 
contributions can enhance the cross section obviously, that is to say cr^^^^ > '^1^^^ i while 
in other regions, we have o'^^^^^ < (^g^^'^ ■ The figure shows that when we have small and 
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moderate tan (3 values, the scalar quark loop contribution to the A^Z^ associated production 
at the LHC is most obvious. We shall also see later from Fig. 7 that when tan/3 has small or 
moderate value, the contributions from gluon fusion subprocess is dominant. Therefore, it is 
possible to use the experimental measurement of the APZ^ associated production at the LHC 
to disentangle the MSSM from the general 2HDM in these parameter space regions. 

In Fig. 5, the cross section of the parent process pp ^ gg ^ + X at the LHC versus 

tan/3 is plotted. The values of the neutral CP-odd Higgs boson A^ mass are set to be 200 GeV, 
400 GeV and 600 GeV, respectively. Prom the figure, we also find that scalar quark contributions 
can either enhance or suppress the cross section of the parent process as shown in Pig.4. Pig.5 
together with Pig.4, show that when the value of niA is greater than 400 GeV and tan /? < 8, the 
contribution from the scalar quark loop diagrams increases with the decrement of tan (5. In Fig. 5 
the two curves for ruA = 400 GeV and 600 GeV demonstrate that when tan/3 < 8, the scalar 
quark contribution suppresses the cross section, which means o-^^^^ < (^1^^^ i ^■nd while the 
scalar quark contribution enhances the cross section when tan/3 > 20. These features can be 
also seen from Fig.4. The curve for tan /3 = 2 in Pig.4, demonstrates that f^^g^^^ is about two 
third of CT^^^^ quantitatively, while the curve for tan/3 = 32 shows cr^g^^^ > CTgg^^- 

Pig.6 displays the differential cross section da/dpT of the process pp ^ gg ^ A^Z^ at the 
LHC versus transverse momentum pT with ^/s = 14 TeV and the pseudo-rapidity being in 
the range of \r]\ < 2. The A^ mass is set to be 350 GeV, and tan/3 is taken as 2, 7 and 32, 
respectively. We find that for tan/3 = 2, the scalar quark contribution suppresses the differential 
cross section da/dpx- For tan/3 = 32 andpT > 100 GeV, the scalar quark contribution enhances 
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the differential cross section. But for tan (3 = 7, the scalar quark contribution can either enhance 
or suppress the differential cross section in different pT regions. 

In Fig.7, Fig.8 and Fig.9, we plot the the cross sections of the process pp ^ bb ^ J^Z^ + X 
and the process pp — >■ J^Z^ + X contributions from both Drell-Yan and gluon-gluon 
fusion subprocesses in the constrained MSSM [a^^ = o^^^^ + cr^^^^), as the functions of 
the CP-odd Higgs boson mass, the ratio of the vacuum expectation values tan/3 and the 
transverse momentum px, respectively. In these three figures, the full-lines are for the cross 
sections or differential cross sections of the process pp ^ bb ^ J^Z^ + X via Drell-Yang 
subprocess, the dotted-lines are for the process pp A^Z^ + X via both Drell-Yang and gluon 
fusion subprocesses. With the comparison between the a^^^ (or da^^^dpr) and a^^^^ (or 
, we can know in which parameter space in the constrained MSSM, the contribution 
from the loop mediated subprocess gg A^Z'^ is important. 

Fig.7 displays the cross sections oipp bb ^ A^Z^+X and pp — *■ A^Z'^+X at proton-proton 
colliders versus the mass of A^ with -s/s = 14 TeV. We choose tan/3=2, 7 and 32, respectively. 
Prom Fig.7, we find that in the region of tan/3 < 7, the contribution of gluon-gluon fusion 
subprocess in the MSSM enhances the cross section, especially when tan j3 = 2 the contribution 
of the gluon-gluon fusion subprocess is about 80% of the total cross section a^'^\ In fact, the 
gluon-gluon fusion subprocess is the most important A^Z^ associated production mechanism in 
this parameter space. From the figure we see also that when tan /3=32 the difference between 
o"^^^) and cr^-^^ is very small, it means that the contribution of the gluon-gluon fusion subprocess 
is negligible in this parameter space. 
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The cross sections of pp ^ bb ^ A^Z'^ + X and pp J^Z'^ + X at the LHC as the functions 
of tan/3 with ^/s = 14 TeV are shown in Fig.8. The mass of Higgs boson is taken as 200 
GeV, 400 GeV and 600 GeV, respectively. Prom this figure we can find also that gluon-gluon 
fusion subprocess enhances the cross section of the J^Z'^ associated production at the LHC, 
and will become a very important production mechanism when tan/3 < 10. In the region of 
tan /3 > 10, the cross section of the J^Z^ associated production at the LHC are in the range of 
1 - 10^ fb. Even the a^^^ can reach 300 fb when tan/3 = 48 and = 200 GeV. So the J^Z^ 
associated production process may be easily observed experimentally if tan /3 is large enough. 

Fig.9 displays the differential cross sections [da jdpT) of pp A^Z^ + X and pp —>■ bb —>■ 
A^Z^ + X at the LHC as the functions of the transverse momentum pT with the pseudo- 
rapidity being in the range of |?7| < 2. We choose ruA = 350 GeV, and take tan /3=2, 7 and 32, 
respectively. Prom this figure we can see that at high pT region, when tan /? < 7, the difference 
between da^^^^dpT and da^^^ /dpr is obvious, even when tan/3=2, the da^^^^^dpr can be less 
than 1% of da^^^ /dpT, which means that the contribution from the pp ^ gg ^ J^Z^ + X 
process is dominant in this parameter space. But when tan/3 = 32, the contribution to the total 
differential cross sections {da^"^ jdpx) is mainly from the Drell-Yan J^Z^ associated production 
subprocess, and the contribution from gluon fusion subprocess is negligible. 

IV Summary 

In this paper, we studied the neutral CP-odd Higgs boson A° production with the association 
of Z^ gauge boson via both Drell-Yan and gluon-gluon fusion subprocesses in the constrained 
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MSSM at the CERN LHC. Numerical analysis of their production rates is carried out with 
some typical parameter sets in the mSUGRA scenario. Our results show that the cross section 
in the MSSM is clearly enhanced by the gluon-gluon fusion subprocess in the parameter space 
with small or moderate tan (3 value, and we should consider the gluon-gluon fusion subprocess 
in this parameter space in the calculation of the A^Z^ associated production at the LHC. We 
compared above results of the process pp ^ gg ^ A^Z^ + X in the MSSM with those in the 
general two-Higgs-doublet model (2HDM), where the cross section of subprocess gg ^ A^Z'^+X 
is contributed only by quark loop diagrams. We find that the contributions from the scalar quark 
loops in the MSSM can either enhance or suppress the cross section obviously and cannot be 
neglected in some parameter space. The results show also that the A'^Z^ associated production 
at the LHC is strongly related to the parameters tan/3 and the mass of A^. The total cross 
section increases with increment of tan/3, and decreases with increment of niA- 
Acknowledgments: This work was supported in part by the National Natural Science Foun- 
dation of China and a grant from the Education Ministry of China . 
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Figure Captions 

Fig.l The relevant Feynman diagrams for the subprocess bb J^Z'^ in the MSSM at the 
tree-level: (a) s-channel diagrams, (b) u- and t-channel diagrams. Note that Fig. 1(b) includes 
the diagram created by exchanging two final states. 

Fig. 2 The relevant Feynman diagrams for the subprocess gg Z^A^ at the one loop-level 
(including only the quark loop diagrams). 
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Fig.3 The relevant Feynman diagrams for the subprocess gg J^Z^ at the one loop-level 
(including only the scalar quark loop diagrams). 

Fig.4 The cross sections o^l^^^ and cr^g^^^ of the process pp ^ gg ^ A^Z'^ + X, as the 
functions of the mass of Higgs boson A^. The input parameter tan/3 is taken as 2, 7 and 32, 
respectively. 

Fig.5 The cross sections cr^^^^ and cr^^^^ of the process pp ^ gg ^ A^Z^ + X, as the 
functions of tan (3. The mass of Higgs boson A^ is taken as 200 GeV, 400 GeV and 600 GeV, 
respectively. 

Fig.6 The differential cross sections dag^^^ /dpT and da^^^^ /dpT of the process pp —>■ 
gg A'^Z'^ + X, as the functions of the transverse momentum px in the mSUGRA scenario 
at the LHC with ^/s = 14 TeV, uia = 350 GeV and the pseudo-rapidity being in the range of 
|j7| < 2. The ratio of the vacuum expectation values tan/3 is taken as 2, 7 and 32, respectively. 

Fig.7 The cross sections a^^^^ and cr^^) of the process pp —>■ J^Z^ -|- X as the functions of 
the mass of Higgs boson The ratio of the vacuum expectation values tan/3 is taken as 2, 7 
and 32, respectively. 

Fig.8 The cross sections a^^^^ and a'^^ of the process pp —>■ A^Z^ -|- X as the functions of 
tan /3. The mass of Higgs boson A° is taken as 200 GeV, 400 GeV and 600 GeV, respectively. 

Fig.9 The differential cross sections da'^^^^ / dpT and da'^^^ /dpT of the process pp — > A^Z*^ + 
X, as the functions of the transverse momentum pT in the mSUGRA scenario with -s/s = 14 TeV, 
ruA = 350 GeV and the pseudo-rapidity being in the range of |r/| < 2. The ratio of the vacuum 
expectation values tan/3 is taken as 2, 7 and 32, respectively. 
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Figure. 1 



